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ABSTRACT 

The halo structure at high Galactic latitudes near both the north and south poles 
is studied using SDSS and SuperCOSMOS data. For the south cap halo, the archive 
of the SuperCOSMOS photographic photometry sky survey is used. The coincident 
source rate between SuperCOSMOS data in Bj band from 16 m 5 to 20 m 5 and SDSS 
data is about 92%, in a common sky area in the south. While that in the Rf band 
is about 85% from 16 m .5 to 19 m .5. Transformed to the SuperCOSMOS system and 
downgraded to the limiting magnitudes of SuperCOSMOS, the star counts in the 
northern Galactic cap from SDSS show up to an 16.9 ±6.3% asymmetric ratio (defined 
as relative fluctuations over the rotational symmetry structure) in the Bj band, and 
up to 13.5 ± 6.7% asymmetric ratio in the Rf band. From SuperCOSMOS Bj and 
Rf bands, the structure of the southern Galactic hemisphere does not show the same 
obvious asymmetric structures as the northern sky does in both the original and 
downgraded SDSS star counts. An axisymmetric halo model with n=2.8 and q=0.7 
can fit the projected number density from SuperCOSMOS fairly well, with an average 
error of about 9.17%. By careful analysis of the difference of star counts between the 
downgraded SDSS northern halo data and SuperCOSMOS southern halo data, it is 
shown that no asymmetry can be detected in the south Galactic cap at the accuracy 
of SuperCOSMOS, and the Virgo overdensity is likely a foreign component in the 
Galactic halo. 
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1 INTRODUCTION 

The modern use of star counts in the study Galatic structure 
began with Bahcall & Soneira (1980). In Bahcall's standard 
model, the structure of the Galaxy is assumed to be an ex- 
ponential disk and a de Vaucouleurs spheriodal halo. A lot 
of work has been done to constrain and examine this the- 
oretical model, as summarised in Xu, Deng & Hu (2006) 
(XDH06 here after), most of them using only a small sky 
area. The global structure of the Galactic halo can only be 
inferred by different observations of small sky areas with 
different magnitude limits, photometric passbands and dif- 
ferent original observational goals. SDSS provides us with 
the opportunity to examine the large scale structure of the 
Galaxy from optical photometry thanks to its deep photom- 
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etry and large sky coverage. From SDSS observational data 
it is clear that the stellar halo of the Galaxy is aymmetric, 
contrary to what has been generally assumed. From colour 
star counts it is obvious that the asymmetric projected stel- 
lar number density is produced by halo stars. There are two 
possible explanations for such a halo structure. Firstly, that 
there are some large scale star streams embedded in the axi- 
symmetric smooth structure of the Galactic halo ( Juric et al. 
2005). Secondly, that the galactic stellar halo is intrinsically 
not axi-symmetric (Newberg & Yanny 2005; XDH 2006). 
Based on the data we have so far, some combination of the 
two might also be possible . In Paper I, we tested the second 
option and fitted the observational data with triaxial halo. 
The triaxial halo model fits fairly well the projected number 
density near the northern cap of the Galactic stellar halo. 
However, in some sky areas, the triaxial halo model cannot 
reproduce the actual star counts. The multi-solutions that 
are intrinsic in fitting the observational data with the triax- 
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ial halo model make the interpretation of the data somewhat 
difficult. On the other hand, the alternative option where the 
asymmetry of the halo is caused by large scale star streams 
also has some problems, even if the overwhelmingly large 
Virgo overdensity that covers nearly a quarter of the north- 
ern hemisphere can be explained by a large scale star stream. 
The observed underdensity near Ursa Major with respect to 
the assumed axi-symmetric halo still challenges such a pic- 
ture. Nevertheless, the conservation of such a huge structure 
in the gravitational well of the Galaxy certainly needs to be 
verified. Martmez-Delgado et al. (2006) show that the Virgo 
overdensity can be reproduced by the dynamical evolution 
of the Sgr stream. Assuming a certain structure of the stel- 
lar halo (an oblate ellipse), their numerical simulation can 
predict an overdensity on a few hundred square degree scale. 

Although it is the most advanced photometric sky sur- 
vey in terms of depth and data quality, SDSS does not have 
good data coverage near the southern cap of the Galaxy 
which is, of course, crucial in understanding the overall 
structure of the stellar halo. Limited to the sky coverage 
of SDSS photometry database, it is probably premature to 
draw a firm conclusion on the stellar halo structure. Assum- 
ing that the Galactic stellar halo is non axi-symmetric, and 
can be described by a triaxial model, there must be some 
corresponding evidence in the southern hemisphere similar 
to what is found in XDH06 for the northern cap. In the ax- 
isymmetric halo model, the maximum star counts should be 
at longitude I = 0° (due to the location of the observer). In 
the case of a triaxial halo, however, the maximum projected 
number density also depends on a certain parameters of the 
halo including azimuth angle, axial ratios and the limiting 
magnitudes of the observations. In the simplest case, the 
plane defined by the primary and the middle axis of the tri- 
axial halo stays in the Galactic disk, the azimuth angle is 
only related with the angle between the primary axis and 
the direction of the Galactic center from the Sun, therefore 
the expected star counts and asymmetric ratio of northern 
and southern sky ought to be mirror symmetric with respect 
to the Galactic plane, i.e. what was found in the north cap 
should also be found in the south under such a halo model. 
If the two planes do not overlap, the situation will be more 
complicated, but similar results should still hold. 

It is also interesting to examine archived sky survey 
data that has the good coverage and reasonable quality in 
the southern Galactic halo: the photographic photometry 
of SuperCOSMOS is ideal for this purpose. As reviewed 
by Hambly et al. (2001a), photographic observations for 
the Galaxy started in the late nineteenth century. In the 
1930s, the development of Schmidt telescopes with wide 
fields of view further advanced photographic surveys. The 
1.2-m Palomar Oschin, 1.0-m ESO and 1.2-m UK Schmidt 
telescopes finished the photographic whole sky survey in the 
last century, such surveys form a legacy library for examin- 
ing the structure of the Galaxy. In the late twentieth century, 
the photographic plates were eventually digitized using mi- 
crodensitometry and digital electronics machines. There are 
several major programs to digitize the photographic plates, 
of which SuperCOSMOS is one. In Hambly et al. (2001a), 
a general overview of these programmes (APM, APS, COS- 
MOS, DSS, PMM, SuperCOSMOS) is presented. The dig- 
itized photographic sky survey of SuperCOSMOS provides 
a catalog of three bands, namely blue(_Bj), red(Rp) and 



near- infrared (Ivn), which have deeper detection limit for 
the same detection completeness compared to other similar 
survey programs (see fig2 of Hambly et al. 2001a). We there- 
fore adopt the SuperCOSMOS data archive for our present 
study. 

In section 2, the observational data are described and 
the stellar source cross identification between the Super- 
COSMOS data and the SDSS data is carried out, and the 
viability of using SuperCOSMOS data to study the struc- 
ture of southern Galactic stellar halo is discussed. In section 
3, downgraded SDSS and SuperCOSMOS observational of 
star counts results are presented. In section 4, the model 
fits to the SuperCOSMOS star counts are introduced. In 
section 5, the SuperCOSMOS observational data and theo- 
retical models are compared, and SuperCOSMOS southern 
sky star counts and SDSS downgraded northern sky star 
counts also compared and analyzed. In section 6, the result 
of star counts is summarized. 



2 THE OBSERVATIONAL DATA 

2.1 SuperCOSMOS photometric data 

The SuperCOSMOS Sky Survey is a digitized photography 
sky survey. It is described in detail in a series of papers by 
Hambly and collaborators (Hambly et al. 2001a; Hambly, 
Irwin & MacGillivray 2001b; Hambly, Davenhall & Irwin 
2001c). 

The SuperCOSMOS photography atlas of the Super- 
COSMOS sky survey includes blue(Bj), red(7?F), and near 
infrared (Ivn) passband photometric survies carried out by 
UK Schmidt survey for -90° < Dec < +2.5°, ESO Red 
Survey of —90° < Dec < —17.5°, and Palomar surveies in- 
cluding, POSS-I Red Survey for -20.5° < Dec < +2.5°, 
POSS-II Blue Survey for -2.5° < Dec < +90.0°, POSS-II 
Red Survey for -2.5° < Dec < +90.0° . Data of Bj band has 
90% or about detection completeness from 16 m .5 to 20 m .5, 
and that of Rf band has same completeness from 16 m .5 to 
19 m .5 mag. The photometric data has a magnitude error of 
m .3, but color is accurate to about m .16(Hambly et al. 
2001a). 

There are two interface applications of SuperCOS- 
MOS: the SuperCOSMOS Sky Survey (http:// www- 
wfau.roe.ac.uk/sss, SSS hereafter) and the SuperCOSMOS 
Sky Archive (http:// surveys. roe. ac.uk/ssa, SSA hereafter). 

Images of small sky areas and catalogs from the Su- 
perCOSMOS sky survey can be downloaded. We thank the 
SuperCOSMOS working group who made all the data avail- 
able to the community. The SSA only includes photomet- 
ric data from UKST and ESO. As made clear by Hambly 
et al. (2001a), although the entire sky is digitized, the data 
in this archive is released progressively. The total amount 
of data is enormous, only F type stars (selected by 0.504 
< Bj - R F < 0.8236 ) from 20 m .4 to 20 m .415 are adopted 
to show the sky coverage which is in the upper panel of fig- 
ure 1. The survey covers most of the high latitude southern 
sky, and a little of the northern hemisphere. The clump at 
(l,b)=(302.616°, -44.580°) is the SMC, and the clump at 
(l,b)=(280.085°, -30.430°) is the LMC. 

Except for the difference in sky coverage, the two in- 
terface applications use different selection standards. The 
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SSA SQL selection is much more configurable (private com- 
munication by email with Hambly) than that of the SSS. 
For example, there are 4 kinds of B magnitude in the SSA, 
namely classMagB (B band magnitude selected by B image 
class), gCorMagB (B band magnitude assuming the object 
is a galaxy), sCorMagB(B band magnitude assuming the 
object is a star), classB (image classification from B band 
detection) . The most appropriate attribute for point sources 
is sCorMag, while the most possible class of an object from 
all three bands is provided by parameter "meanclass". The 
SSS only includes selection parameters applied to the pri- 
mary passband, coresponding to classB of the SSA in the 
example. Our aim is to count the stars in each selected sky 
area, and using classB will lose some stars due to not synthe- 
sizing information of all the three bands. This will influence 
the result of star counts seriously. So "sCorMagB" of the 
SSA data is selected to carry out the study and the "mean- 
class" is limited to 2(star label). Because the SSA only covers 
limited sky areas of high latitude northern galactic hemi- 
sphere (upper panel of figure 1) we cannot directly compare 
SuperCOSMOS star counts of northern sky with those of 
SDSS. 

The SSA includes R F band data from both UKST and 
ESO. However data from the Rf band of UKST is deeper 
than that of ESO. Therefore, only UKST is adopted. The 
detailed instrumental specifications of UKST can be found 
in Cannon (1984), the main parameters of the survey tele- 
scope and instruments are listed in table 1. 

As demonstrated in the upper panel of figure 1, the 
UKST atlas of SuperCOSMOS covers most of the high 
Galactic latitude southern hemisphere. The structure of the 
Galactic halo near the southern cap can be studied using a 
stellar photometry catalog selected in a similar way as we 
did for the northern sky in XDH06, shown here in the lower 
panel of Fig 1. The selected sky area for this work is shown 
in lower panel of figure 1, the Lambert projection of south- 
ern hemisphere. Each of the squares represents a rectangular 
sky area of about 2° x 2° . Some of the selected sky areas may 
be trimmed if sitting near the survey's edge, or the region 
is masked by contaminats such as saturated bright stars, 
or clumps such as the dwarf galaxy IC1613 in (130°, —60°). 
The first group of sky areas are along a circle of b — —60°, 
equally spaced by 10°. The other 12 groups are a selection 
of sky areas along longitudinal directions equally spaced by 
30°. At a given longitude, the sky areas are selected by a 
step of 5°. This selection of sky areas can evenly cover the 
southern Galactic cap, so that the global structure of the 
halo near the southern pole can be examined. 

2.2 Cross checking of SuperCOSMOS and SDSS 
data sets 

In our previous work (XDH06), SDSS data is used to 
study the structure of the Galactic stellar halo near 
the Northern Galactic pole, The SDSS catalog provid- 
ing a uniform and accurate photometric data set. The 
five broadband filters, u,g,r,i,z are 95% complete to 
22 m .0,22 m .2,22 m .2,21 m .3,20 m .5 respectively, and the un- 
certainty in the photometry is about 3% at g=19 m (Chen et 
al. 2001). 

Compared to the high-quality photometry data of 
SDSS, the SuperCOSMOS data has a narrower dynamic 



range, lower magnitude limit and larger photometric error, 
due to photographic photometry. To evaluate any uncertain- 
ties due to misclassifications and the relatively less accurate 
photometry of SuperCOMOS, a comparison in areas com- 
mon to both surveys is needed. 

The photometric calibration between SDSS and Super- 
COSMOS has been made available by the 2dF Galaxy Red- 
shift Survey(2DFGRS) Final Data Release Photometric Cal- 
ibration which defines a set of color equations in its final 
dataQ. The Bj band is correlated with SDSS g and r band, 
Bj = 0.15 + 0.13 X (g-r), while Rf is very similar to SDSS 
r band, 7?F=r-0.13. The results of such color calibration are 
shown in figure 2 . 

The two small sky areas with superpositions of SDSS 
and SuperCOSMOS surveys in both the northern and south- 
ern sky are chosen to examine the color equations and 
the classification of SuperCOSMOS objects. The north- 
ern area is located at (/, fe)=(280° , 60°) with 2° x 2° field 
of view (FOV), the southern area is at {1, 6)=(62°, -59°) 
with 1° x 4° FOV. The equinox of SuperCOSMOS data 
associated with the photometric image library is J2000.0. 
The position accuracy of SuperCOSMOS is ±0.2arcsec at 
Bj = 19 m , R F = 18 m , ±0.3arcsec at Bj = 22 m ,R F = 21 m 
(Hambly et al 2001c). Taking into account proper motion, 
cross-identification is carried out between SDSS and Super- 
COSMOS in the two superimposed sky areas in a identifica- 
tion criterion box of m .3 and 10 arcsec. In such a box, mul- 
tiple sources can be present, the pair of stars with the near- 
est coordinates and magnitudes are identified as the same 
source. We take the SDSS data as the "true" values of both 
position and magnitude. Based on the matched star list in 
the two areas, uncertainties in the magnitude of SuperCOS- 
MOS photometry for each object can be measured. The sys- 
tematic error calculated this way infers the error of the color 
equations from 2DFGRS calibration; while the scatter can 
be used to measure the error in SuperCOSMOS photometric 
data. Fitting the systematic error with a 2nd order polyno- 
mial, the color equations are refined. Using the modified 
color equations, Bj and Rf magnitude of SDSS data is de- 
fined as BjsDSS=g+0.15+0.13x(g-r)+ Amod, Rf.sdss=?- 
0.13+Amod. Iterating the cross-identification procedure re- 
duces the systematic error. The error in the SuperCOSMOS 
data in the Bj band is found to be eg, = B.j - Bjsdss, 
and that in Rf £Rj = Rf - Rfsdss- The variance of the 
errors as functions of magnitude is obtained from fitting the 
scatter with a gaussian. 

After such modification, and repeating the cross- 
identification, the source matching ratios between the two 
surveys are improved. In the end, the SuperCOSMOS data 
matches that of SDSS in the Bj band magnitude limits by 
92-93% in a 10 arcsec and m .3 box. For the R F band, 
the matching ratio can be raised to 85% or larger in the 
16 m .5 — 19 m .5 magnitude range. The matching ratio of the 
Rf band is not as good as that of the Bj band, this is 
likely due to the lower sensitivity in the Rf band. A 85% is 
still lower than the intrinsic completeness estimated for dif- 
ferent surveys in the SuperCOSMOS atlas (see figure 12b of 
Hambly et al. 2001b). This is possibly caused by the brighter 
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Table 1. Parameters of telescope of UKST 



site Siding Spring Mountain, -31°S 

aperture 1.24m 

focal, focal ratio 3.07m, f/2.5 

photographic plates Kodak Illa-J emulsion, 356mm 2 , 67.1 arcsec mm -1 , 6.5° x6. 5° 

primary pointing accuracy ±6 arcsec r.m.s. 



magnitude limit of the SSA compared to that of SDSS as the 
bright stars are saturated, therefore infiuencin more neigh- 
bours. 

In the upper and lower panel of figures 3, the contours in 
the color-color diagram represent the SuperCOSMOS data 
in the three bands that are cross-identified in the SDSS data. 
Black points over-plotted on the contours are the matched 
stellar sources, while the crosses represent SuperCOSMOS 
sources which are unmatched. 



3 OBSERVATIONAL STAR COUNTS 

3.1 Star counts from downgraded SDSS data 

The examination of halo structure through star counts de- 
pends critically on the depth of the photometry. The SSA 
has a narrower dynamic range and shallower detection limit 
than SDSS. A test is carried out to check if the asymmet- 
ric structure found in XDH06 is still present with the shal- 
lower limit of SSA data. The data used in XDH06 is down- 
graded by applying the SuperCOSMOS magnitude limits, 
photometric errors of SuperCOSMOS are also added to the 
SDSS data. A Monte-Carlo method is used to reproduce the 
photometric errors as of SuperCOSMOS en , , en F (Rock- 
osi private communication). Gaussian errors similar in size 
to those of the SuperCOSMOS data are added to the mag- 
nitude of each star before measuring the star counts. We 
find that the results of XDN06 are recovered, with the av- 
erage fluctuation raised only by about 3.7%. After trans- 
forming into the SuperCOSMOS system, the errors are 
16 m 5 < Bjsdss < 20 m .5 and 16 m .5 < Rfsdss < 19 m .5 
respectively. 

Figure 4 and Figure 5 show the star counts from the 
SDSS data with same sky areas as in XDH06 but down- 
graded to the SuperCOSMOS magnitude limits, for Bjsdss 
and Rfsdss respectively. From the present SDSS public 
data release, the sky area I = 210° is now added. Panel a) 
is for star counts in sky areas along the b = 60° circle. Pan- 
els b)-f) are for star counts of sky areas along the selected 
longitudinal directions paired by mirror symmetry on the 
both sides of the I = 0° meridian. The asymmetric struc- 
ture still appears clearly with the magnitude limit of the 
downgraded SDSS data (especially in figure 4a). The asym- 
metric structure is not so prominent as with the original 
SDSS magnitude limits, but we can still see that the star 
counts in / 6 [180°, 360°] are systematically higher than in 
I £ [0°, 180°]. The largest asymmetry of star counts appear 
in panels b), c) and d). In panels e) and f), the errors are so 
large at the downgraded limits that the asymmetric differ- 
ences between sky areas found in XDH06 are only marginally 
visible. As in figure 4, figure 5 shows the results of star 
counts from the Rfsdss data. Again, the most prominent 
excess over mirror symmetry is found in panels b) , c) and d) . 



However, the Rfsdss band magnitude limit is fainter than 
that of the Bjsdss band, which leads to weaker features of 
asymmetry than figure 4. Tables 2 and 3 describe the asym- 
metric ratio and its uncertainty in the downgraded SDSS 
data. Columns 1-4 are the Galactic coordinates (I and 6), 
counted numbers and the corresponding errors for sky areas 
with I 180°, and columns 5-8 are the same quantities for 
sky areas on the other side of the I — 0° meridian. Com- 
parison is between sky areas paired with mirror symmetry 
with respect to the I = 0° meridian. The asymmetric ra- 
tios are defined by: asymmetricratio = (number density-?, — 
number density i) / (number density i) x 100% which are given 
in column 9; column 10 gives the uncertainties in the ra- 
tios which are inferred from the error of the number den- 
sities (tables [4]-[5] all have the same entries, but for dif- 
ferent data). The asymmetric ratios measured from down- 
graded SDSS data are all positive with one exception which 
is very near to zero, this means that all the sky areas in 
/ 6 [180°, 360°] have higher projected number densities 
than those in I £ [0°,180°]. The largest asymmetric ratio 
is 16.9 ± 6.3% in the Bjsdss band and 13.5 ± 6.7% in the 
Rfsdss band. 

Therefore if there are similar levels of asymmetric struc- 
ture in the southern Sky, they should be visible even with 
the SuperCOSMOS magnitude limit. 

3.2 Southern Galactic cap: star counts from 
SuperCOSMOS data 

In XDH06, star counts from SDSS data show a prominent 
asymmetric structure in the northern Galactic hemisphere 
through comparing the projected number densities of sky 
area pairs with mirror symmetry on both sides of the I — 
0° meridian. We will use the same method to examine the 
structure of stellar halo in the southern sky, in particular to 
check whether the halo structure has the same features or is 
different from its northern counterpart. 

Star counts for southern sky from the SuperCOSMOS 
data are shown in figures 6 and 7. Star counts in each sky 
area are plotted using triangles and squares. Panel a) shows 
the results of star counts for the selected sky areas along a 
circle of b = —60°. Panels b), c), d), e) and f) are for sky 
areas along the longitudinal directions, also paired with mir- 
ror symmetry on the both side of the I = 0° meridian. Each 
of the sky areas is divided into four subfields to account for 
the fluctuation of star counts over the average value of the 
area. The fluctuations calculated for all sky areas this way 
are used as error bars in the plots. The error bars actually 
measure the intrinsic fluctuations of the projected number 
density, and the uncertainties in classification and photom- 
etry. The average error of star counts will be discussed in 
section 6. 

Dividing the southern cap into two halves by the I = 
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Table 2. The relative deviations of the sky field pairs for Bjsdss£ [16 m .5, 20 m .5] mag. 
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0°, 180° meridian, the data for both Bj and Rf bands show 
that the structures of the two halves are basically symmetric 
within statistical errors. This is clearly shown in panel b)-f) 
of figures 6 and 7 . 

The Bj band data shows smaller error bars and obvious 
smoother structure in the projected number density distri- 
bution than the Rf band data does. Sizable fluctuations 
over a axisymmetric structure do exist in the Bj band data. 
In two pairs of data, i.e. (150° , -60°)& (210°, -60°) and 
(90°, -60°)&(270°, -60°), the projected number density at 
I > 180° side is higher than the other side (/ < 180°). While 
the pair of (60°, -70°)& (300°, -70°) shows a reversed ex- 
cess. 

The star counts from the Rf band have a larger scatter 
than that of the Bj band. The Rf band data also has less 
coincidence in classification with SDSS data than the Bj 
band data. Moreover, its limiting magnitude is shallower 
than the Bj band by about one magnitude. For example, 
for FO-type stars, the distance limits given by the Bj band 
are from 5.23Kpc to 32.98Kpc, while those defined by the 
Rf band are from 6.46Kpc to 25.72Kpc; and for F8-type 
stars, they are 2.68Kpc to 16.89Kpc for the Bj band, and 
3.79Kpc to 15.09Kpc for the Rf band. Selecting redder stars 
from a shallower box in the Rf band, star counts show larger 
deviations. 

In both the Bj and Rf bands, there is a odd data 
point at (130°, —60°), which has a projected number density 
obviously lower than its neighbor sky areas. Because this sky 



area is near the edge of the survey, it is very likely that this 
is a boundary effect. 

Table 4 (for Bj) and table 5 (for Rf) list the pro- 
jected number densities and their corresponding errors, and 
the asymmetric ratio measured in the SuperCOSMOS data. 
Comparing the uniform positive asymmetric ratios in the 
downgraded SDSS data (tables 2 and 3), the values given 
by SuperCOSMOS are quite irregular, with apparantly ran- 
dom positive and negative values. 



4 THE THEORETICAL MODEL 

From SuperCOSMOS star counts there is no obvious asym- 
metric structure in the southern hao. A theoretical axisym- 
metric halo model is therefore adopted here. Because the Rf 
band data is shallower and less consistent with SDSS data, 
only Bj band data is used to constrain the model parame- 
ters. 

From the Bahcall standard model (Bahcall & Soneira 
1980), the projected number density in a certain apparent 
magnitude interval along a fixed direction can be described 
by the integral of a density profile and a luminosity function. 

A(m 1 ,m 2 ,e,b) = / dm / R 2 ARp{v) 4>{M) dft, (1) 

J mi JO 

where mi , 7712 are the limits of the given magnitude in- 
terval, R is the heliocentric distance of a star, p is the density 
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Table 3. The relative deviations of the sky field pairs for Rfsdss& [16™. 5, f9 m .5] mag. 
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profile of each stellar population, and 4>(M) is the luminosity 
function of the population. 

For the thin and thick disk components, the density 
profile is assumed to be exponential, 

p(r) =exp[-\z\/H -x/h], (2) 

where \z\ is absolute value of height of a star above the 
Galactic plane, x is distance between the Galactic centre and 
the projected point of that star on the Galactic plane, H is 
the scale hight, and h is the scale length of the exponential 
disk. 

The power-law halo density profile in Reid(1993) is 
adopted, 

n I n 

PW al + r n ' ' 

where ro is the distance from the Sun to the Galac- 
tic centre, and ao = 1000 is a normalisation constant. 
r is the distance from the star to the Galactic centre. 
r = \Jx 2 + (y/p) 2 + (z/q) 2 as in XDH06. x, y, z define the 
position vector in three axes of coordinate frame adopted 
in paper I. p, q are the axial ratios of the middle axis and 
shortest axis to the major axis respectively. The triaxial halo 
model naturally degenerates to asymmetric halo model when 
p=l The luminosity functions of the halo, thick and thin 
disk components in Bj, Rf bands are transformed from the 
luminosity functions of Robin & Creze(1986), with Bj=B- 
0.304x(B-V) and R F = R + 0.163 x (V - R) as provided 



by the photometric calibration of the final data release of 
2DFGRS. Given the Rf luminosity function, star counts in 
the R F band can also be obtained, but these are not used 
for model fitting for the reasons given above. 

The three dimensional extinction model of the Milky 
Way derived from COBE observations is adopted four our 
model. Directly correcting the observational data for extinc- 
tion is not possible due to the lack of distance informa- 
tion for individual stars (Drimmel, Cabrera-Lavers & Lopez- 
Corredoira 2003), however we solve this problem by applying 
COBE extinction data to the theoretical model. 

Using equation (1), the projected number density of 
each sky area can be obtained. To reveal the distribution 
of star counts in apparent magnitude, the Bj band mag- 
nitude is divided into 8 bins (16 m .5 to 20 m .5, in steps of 
0.5 m ), the number density in each magnitude bin is then 
calculated. Constraining star counts in Bj band magnitude 
limits of 16 m .5 to 20 m .5, a non-neligible number of stars 
have no corresponding Rf band data, therefore color counts 
need to be treated with special care, and that will be dis- 
cussed later in section 5.2. 

As a continuation of XDH06, the present work is focused 
on halo structure near the southern cap of the Galaxy. For 
a better comparison between the present work and XDH06, 
the parameters of the thin and thick disks are fixed with 
the values used in XDH06, which were taken from Chen 
et al.(2001). Only the halo parameters are adjusted to fit 
the SuperCOSMOS observations, the scope of parameters is 
listed in table 5. In an axisymmetric halo model there are 
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Table 4. The relative deviations of the sky field pairs for Bj£ [16 m .5, 20 m .5] mag. 
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only two parameters: n is power law index of halo density 
profile, and q is the axial ratio z/x. 

A x 2 minimization is adopted to compare the theoreti- 
cal results and the observational data sets (Press et al. 1992). 
For a non-Gaussian distribution of discrete data, Pearson's 
X 2 is used, 



X 2 = V (% - g ' )2 (4) 

The meanings of all the symbols are described in 
XDH06. x' 2 is calculated in order to evaluate the similarity 
between the theoretical projected number density and the 
observational data. Xbin describes the difference between the 
distribution of the theoretical star counts in apparent mag- 
nitude bins and that of observations for each sky area. Xu„ 
is the average value of xlin m all sky areas. 



5 RESULTS AND ANALYSIS 

5.1 Fitting SuperCOSMOS star counts with the 
axisymmetric model 

The theoretical projected surface number densities are cal- 
culated using the axisymmetric model described in section 
5, with extinction included. The theoretical model that best 
fits the SuperCOSMOS observational data in both Bj band 
and Rf band is shown in figure 6 and 7 as the solid and 
dashed-lines respectively . The model parameters are n=2.8, 
q=0.7. This is one of the best fitting models in the provided 
parameter space. An axisymmetric model can fit SuperCOS- 
MOS data reasonably well within the statistical error bars. 
The best fit theoretical model (solid line) and the observa- 
tional data (diamonds with error bars)for l=-60° fields are 
shown in figure 6a in which the data shows an irregular pat- 
tern of deviations from the symmetric model. The projected 
number densities at (40°, -60°), (60°, -60°), (270°, -60°), 
(290°, —60°) are higher than the model, while the observa- 
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Table 5. The relative deviations of the sky field pairs for Rp£ [16 m .5, f9 m .5] mag. 
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tional data at (110°, -60°), (340°, -60°), (350°, -60°) are 
lower than the model. It is clear from figures 6b-6f that the 
two theoretical lines do not overlap perfectly due to differ- 
ent extinctions. In figure 6e of b=-70°, 1=120°, the value 
of the theoretical curve has a dip because that sky area 
(120°,-70°) has a large extinction from COBE, such that 
when the distance is larger than 550PC, the extinction takes 
Av — m .174. While in other areas around this, the extinc- 
tion ranges only m .07 - m .08. 

The above model parameter set is used to calculate the 
theoretical Rf band star count, and the results are plot- 
ted in figure 7. In figure 7a, star counts from (40°,— 60°) 
to (110°, -60°), (270°, -60°) to (340°, -60°) also fluctuate 
around the theoretical value. In figures 7b-7f, the pairs of 
(120°, -75°) & (240°, -75°) and (150°, -75°)&(210°, -75°) 
show counts higher than the theoretical line, while counts in 
all other areas are random around the theoretical prediction. 

As well as calculating the stellar projected number den- 
sity, we also compare the theoretical and the observational 



star counts in apparent magnitude bin for each sky area. 
Figure 8 shows, as an example, the observational (gray di- 
amond) and the theoretical (dark line) star counts in 12 
sky areas of b = —60°, the Galactic coordinates of each sky 
area are indicated in the corresponding panel. As shown in 
these plots, the theoretical model can fit observation data 
fairly well, but with a few exceptions. In (90°,— 60°) and 
(330°, -60°) at the bin of 20 m to 20 m .5, the theoretical 
value is higher than the observational one. Similar to what 
is shown in figures 6 and 7, the distribution of star counts 
in apparent magnitude also fits fairly well a homogeneous 
axisymmetric structure. 

Using equation (4), \ 2 and Xu n f° r each parameter grid 
can be obtained. The contour plots of \ 2 and xlin m tn e n-q 
plane are presented in figure 9. The minimum value of \ 2 
is 1.53, and the maximum is 11.915, while the values for 
xl in are 0.939 and 3.424 respectively, 20 levels of contours 
are used. The inner-most (smallest values of \ 2 and xti„) 
contour indicates the best fitting combinations n and q. The 
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Table 6. Input range of parameters of theoretical model 



Parameter 


Lower limit 


Upper limit 


step 


11 


2 


4 


0.1 


q 


0.4 


1.0 


0.1 



Table 7. The best-fiting parameters 
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q 




x 2 
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.714 


0.990 


2 


.6 


0.6 
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.658 
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0.948 


2 


.8 
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open diamonds in both panels of figure 9 indicate the most 
favorable parameters given by both \ 2 an d xL minimiza- 
tions, which are listed in table 7. The fitting of the observed 
projected number density using one of the best combinations 
(n=2.8, q=0.7) is shown in figure 6. 

5.2 Comparison between star counts of the 
northern and southern Galactic caps 

In the previous subsection, the southern sky projected sur- 
face number density of SuperCOSMOS Bj band data is fit- 
ted by an axisymmetric stellar halo model. As discussed 
above, star counts of the northern sky show asymmetric 
structure due to an excess of halo stars for I > 180° (see 
XDH06 for details). The presence of the same feature in the 
southern sky is the main concern of this paper. 

To answer this question, we need to compare the distri- 
bution of number density in the north from the downgraded 
SDSS data and that of SuperCOSMOS data in the south. 
Singling-out the halo population from star counts is now 
required. With the data we have, the halo and disc pop- 
ulations can only be roughly distinguished through colors 
based on photometric data. SuperCOSMOS Rf band data 
has only an 85% coincidence with SDSS data, which makes 
our analysis somewhat less accurate. However, this factor 
only affects the total number of stars that can be used in 
statistics in color, and will raise the level of random error in 
the final result. Further to this aim, Rf band data is still 
again used to obtain the star counts in color. 

Figure [10] shows the projected number density of SDSS 
downgraded data of 6 = 60° and SuperCOSMOS data of 
b — —60°. Both data sets are constrained by Bj and Rf 
band magnitude limits (16 m .5 < B JS dss,Bj < 20 m .5, 
16 m .5 < Rfsdss,Rf < 19 m .5). Black points and gray 



points represent SuperCOSMOS data for b = —60° and 
SDSS downgraded data for b = 60° respectively. To show 
the difference clearly, a 6th order polynomial function is 
used to fit for each data set. The SDSS downgraded data are 
systematically higher than SuperCOSMOS data. There are 
two possible reasons for this: firstly, a systematic deviation 
between the two systems; secondly, an intrinsic difference 
between the north and the south. From I = 0° to 240°, the 
two curves have similar shape, showing a possible system- 
atic deviation between the two systems. While from / = 240° 
to 360°, data set for b — 60° shows an obvious excess over 
that of b = —60° after considering the systematic deviation. 
The largest excess appear around I = 330°, coincident with 
the Virgo overdensity (Newberg & Yanny. 2005; Juric et al. 
2005; XDH06). 

Figure [TT] shows the projected surface number density 
in Bj — Rf color space for (90°, 60°) (the gray line of up- 
per panel) & (270°, 60°) (the black line of upper panel), and 
(90°, 60°) (the gray line of lower panel) & (90°, -60°) (the 
black line of lower panel). In the upper panel, the distri- 
bution of SDSS downgraded data in color shows the same 
property as that in XDH06, the halo populations (blue peak) 
in the sky areas I > 180° have an excess over those I < 180° , 
while the disk populations (the red peak) are basically the 
same. The Lower panel shows that both the SDSS down- 
graded data and the SuperCOSMOS data sitting at two op- 
posite sides of the Galactic plane have a double peak struc- 
ture in color space. The disk population in the two sky ar- 
eas has similar number density while the northern sky star 
counts of halo population have larger numbers than those 
in the southern sky. In figure 1101 the systematic deviation 
between the two curves is caused by the difference in photo- 
metric sensitivity limits between the two systems. The rea- 
son is quite straightforward: the fainter stars between the 
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photometric limits of SuperCOMOS and that of SDSS are 
surely absent from SuperCOSMOS statistics, while possibly 
being present in SDSS catalog. 

The lowest number density of star counts in color ap- 
pears for Bj — Rf = 1.6. The disk population and the 
halo population can be roughly separated by this limit ( 
Bj — Rf > 1.6 for the disk population, and B.j — Rf < 1.6 
for the halo population). Figure [12] demonstrates the differ- 
ence between the selected populations in the north (down- 
graded SDSS data) and in the south (SuperCOSMOS data). 
The upper panel of figure [12] shows the difference between 
the density of the halo population of sky areas along the 
b — 60° circle and that of the b — —60° circle. The lower 
panel is the same as the upper one but for the disk popu- 
lation. It is clear that the difference in disk population in 
the lower panel (the southern Galactic cap) has a random 
distribution around 0, the amplitude of such fluctuations is 
lower than about 40 with no systematic feature; while the 
difference in halo population in the upper panel has obvi- 
ous features of over 200. The systematic deviations between 
the SDSS downgraded data and the SuperCOSMOS data 
are clearly caused by halo stars, i.e. the halo population in 
the north has a certain amount of excess over that in the 
south. Clearly, there is a prominent excess in the range of 
I = 300° -360°. This shows that there is an overdensity only 
in the north, while no such features are found in the southern 
SuperCOSMOS data. 



6 DISCUSSION AND CONCLUSIONS 

From SDSS data covering the northern cap, it has been 
found that the northern halo is not axisymmetric (Newberg 
& Yanny 2005; Juric et al. 2005; XDH06). This feature is 
also visible at shallower magnitude limits (ie. closer halo 
stars) in SDSS data downgraded to the limit of SuperCOS- 
MOS. The main goal of this work was to examine the halo 
structure near the southern cap of the Galaxy using Super- 
COSMOS data. We show that the southern halo structure 
does not have a similar asymmetry to the northern galactic 
cap for the same magnitude limits. 

In XDH06, using very deep SDSS photometry from 
15mag- 22mag, the asymmetry ratio goes up to 23%. The 
magnitude limit of SuperCOSMOS data is from 16 m .5-20'" .5 
for Bj band and 16'". 5-19'". 5 for R F band. Converting the 
SDSS data to the same photometry system and consider- 
ing Bjsdss in the same magnitude range, the asymmetric 
structure is weakened but still detectible, as demonstrated 
by the asymmetry ratios and their errors in tables [2] (also 
see [3] for Rfsdss), the asymmetric ratio only pickes up to 
16.9% ±6.3% . From SuperCOSMOS data in the south, star 
counts shows no asymmetry feature, as shown in tables [4] 
(also [5] for Rf), this is of course linked to the uncertain- 
ties in the data. The RMS is over 7.8% for asymmetry ratio 
measured in Bj. 

Concerning the error of star counts of SuperCOSMOS 
Bj band data, there are three sources contributing. Firstly, 
the SuperCOSMOS data of Bj band has 92~93% identifica- 
tion rate when cross-correlating with SDSS data, this gives 
a error of 8%, at the worst case, in number counts. Sec- 
ondly, the SuperCOSMOS data has an overall photometry 
uncertainty of m .3, which creates an error of 3.7% in the 



results of star counts, as derived from Monte Carlo simula- 
tions. Thirdly, the SDSS photometry is far more accurate 
than that of SuperCOSMOS, therefore it can be regarded 
as the precise system to compare with the later one. There- 
fore we assume that the statistical fluctuations measured in 
XDH06 are the intrinsic stellar density fluctuations in the 
halo, which are 2.53% on average for number counts. Putting 
these factors together, we can estimate the average error in 
SuperCOSMOS star counts as, 

a = x/8 2 + 3.7 2 + 2.53 2 = 9.17 (5) 

Having such an uncertainty in star counts for Super- 
COSMOS data, and considering the level of asymmetry of 
16.9% ±6.3%, it is not possible to draw a firm conclusion for 
the symmetry issue for the stellar halo near the south cap, 
when there is only SuperCOSMOS data available. 

However, when analyzing the population statistics using 
colors, distinct properties of stellar halo structures in the 
north and south can be found. As shown in figure 1121 the 
halo population shows an apparent excess around I — 330° 
in the north (the upper panel) as from the downgraded SDSS 
data, while the same plot for the south gives only random 
fluctuations of the same level as statistical errors. 

We attempt to fit triaxial halo models to both down- 
graded SDSS and SuperCOSMOS data. By directly applying 
models in XDH06, no good fit can be derived, because no 
obvious overdensity such as the Virgo one in the north is 
found in the south. However, this does not exclude the pos- 
sibility to have a triaxial halo after removing the large scale 
star streams. Due to large photometric uncertainties and low 
sensitvity of SuperCOSMOS, an error in star counts around 
9.17% prevents us from making a clear conclusion on this 
point. 

Therefore, the present work can be concluded as the 
following: 

(i) SuperCOSMOS data (SSA) has been used to study 
the structure of stellar halo covering the southern Galactic 
cap. Direct star counts reveal that the structure can be fitted 
by a axi-symmetric halo model. Limited by the photometric 
error and depth of the survey, no asymmetry can be detected 
by star counts. 

(ii) An Asymmetric structure, very similar to what have 
been found using SDSS survey data (Newberg & Yanny 
2005; Juric et al. 2005; XDH06) can be detected by down- 
grading SDSS data to the limiting magnitudes and photo- 
metric error of SuperCOSMOS. 

(iii) A halo population excess, defined by (Bj — Rf < 
1.6), is responsible for the asymmetry structure found in 
the north in downgraded SDSS data, as revealed by both 
direct star counts (figures (4H5] I10[) and statistics in color 
(figures [TT] and I12|) . While for the southern cap, no such 
features are present. 

(iv) Considering the overall symmetry of the Galactic 
halo, the asymmetry discovered in the north (the Virgo over- 
density) is likely to be a foreign component in the stellar halo 
of the Galaxy. However, due to a lack of good photometric 
data, an asymmetry in the stellar halo near the south cap 
beyond SuperCOSMOS limits cannot be ruled out. It is still 
an open question if we have a triaxial halo with large scale 
star streams embedded. 
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(v) For the structure of stellar halo near the southern cap, 
SuperCOMOS data cannot go any further. Better quality 
survey data of the SDSS quality is needed to adress these 
issues. 
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Figure 1. Upper panel: the sky coverage of the SupcrCOSMOS archive in Galactic coordinates as shown by F type stars in the Bj 
band from 20 m .4 to 20™. 415. Lower panel: Lambert projection of the sky coverage of SuperCOSMOS (shown by tiny dots of the same 
selection of stars in the upper panel) , and the sky areas selected for this study (squares) . 
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Figure 2. The color transformation between SuperCOSMOS Bj, Rp bands and SDSS g,r bands. The obvious offset between the two 
systems infers that a systematic correction is needed. See text for details. 
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Figure 3. The cross identification between the SuperCOSMOS and SDSS data sets in color-color space. The contour and the black 
points on the color-color diagram represent the Bj band matched sources in the same overlapping sky areas, limited in m .3 in magnitude 
and 10 arcsec in angular distance box and a magnitude interval of 16 m .5 < Bj, Bjsdss < 20 m .5. The crosses represent the unmatched 
sources. The upper panel: The northern overlapped sky area around (280°, 60°); The lower panel: the southern one around (62°, —59°). 
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Figure 4. Projected surface number density of SDSS data in the same selected sky areas as in XDH06, but downgraded to 16 m .5 < 
BjSDSS < 20 m .5. Panel a), is for the selected areas along a circle of b = 60°, the horizontal axis is the Galactic longitude in degrees; 
while the others arc for the ones along different paired longitudinal directions, with the longitudes indicated in the inlet of each panel, 
all the horizontal axis in panels b)-f) are the Galactic latitude in degrees. 
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Figure 5. The same as figure [4] but for the projected surface number density of SDSS data in Rjsdss downgraded to 16" 
Rfsdss < 19 m -5 
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Figure 6. The same as figure|4] but for the fitting of the surface number density counted from SuperCOSMOS Bj band data. The solid 
and dashed lines are the theoretical predictions, while the diamonds and triangles show the observational data. 



18 Y. Xu, L. C. Deng and J. Y. Hu 



2000 
1500 
1000 

500 


( 

1800 
1500 
1200 
900 

600 
300 





















b=-60 


i a) 





60 120 180 240 300 360 
I 

















» 1=60 




a 1-300 




2500 



g 2000 



"D 

o 

Q_ 



"D 

Q_ 



in 
c 
<D 
"D 



E 

C 

"O 

<D 

O 

'cT 
it 



1500 




1000 



1600 
1400 
1200 
1000 
800 

600 
400 



_ K % * 






' d) 


< 

i 


> 1=90 
1 1=270 




Figure 7. The same as figure [6] but for band results. 
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Figure 8. Model fitting to Bj star counts in a selected sky areas whose Galactic coordinates 
are observational star counts, and the solid lines are the theoretical predictions. 



are indicated in each panel. The grey dots 



20 Y. Xu, L. C. Deng and J. Y. Hu 




power low Index 




power low Index 



Figure 9. The upper panel: Contours of x 2 °f the theoretical models in power law index (horizontal) and axial ratio (vertical) plane. 
The lower panel: Contours of X 2 in of the models in the same plane. The overlapped grids of the minimum contour level are labeled by 
diamonds, which define the best fitting model parameters. 
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Figure 10. A direct comparison between the distributions of the surface number densities of b = 60° sky areas and the b = —60° 
ones, with 16 m .5 < B JSD SS, Bj < 20 m .5, 16 m .5 < R FS DSS, Rf < 19 m .5. The black points and line represent SuperCOSMOS data of 
b = —60° and a polynomial fitting curve. Gray points and line are the corresponding SDSS ones. 
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Figure 11. The upper panel: The color distribution of stars in (l,b) = (90°, 60°) (the gray histogram) and that in (l,b) = (270°, 60°) 
(the black histogram) for SDSS downgraded data with 16 m .5 < Bjsdss < 20 m .5, 16 m .5 < Rfsdss < 19 m .5. An overdensity due to 
halo stars at (l,b) = (270°, 60°) is shown. The lower panel: The same as the upper panel, but for (l,b) = (90°, 60°) (the gray histogram) 
of the downgraded SDSS data and (I, b) = (90°,— 60°) of SuperCOSMOS (the black histogram). An overdensity also due to halo stars 
in the north ((I, b) = (90°, 60°)) compared to its symmetric field in the south ((I, b) = (90°, —60°)) is clear visible. 
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Figure 12. The difference of projected surface number densities between the sky areas along the circles in the north (b = 60°) and south 
(b = -60°). Stars are selected with 16 m .5 < B JS DSS,Bj < 20 m .5, 16 m .5 < Rfsdss,Rf < 19 m .5. The upper panel: Data points arc 
constrained by < Bjsdss ~ Rfsd,SS> Bj — Rf < 1-6 which roughly represents the halo population. The lower panel: Data points are 
constrained with 1.6 < Bjsdss ~ Rfsdss, Bj — Rp < 3.0 which roughly represents the disk population. 



